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Abstract. We calculate transverse-momentum (p t ) spectra and elliptic flow 
of J/^'s in Au-Au and Cu-Cu collisions at RHIC. We employ an earlier con- 
structed 2-component approach for direct and regenerated charmonia which is 
fairly consistent with the centrality dependence of inclusive J/ijj yields at SPS 
and RHIC. Discrepancies with high-pt data in Cu-Cu collisions are addressed 
by implementing effects of a finite formation time and an estimate of B meson 
feeddown, which largely resolve the problem. The sensitivity of the spectra to 
the Cronin effect and different in-medium dissociation rates is studied. The 
predicted elliptic flow in semicentral Au-Au collisions is rather small, ~2-3%. 
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O ■ 1. Introduction 



An "anomalous" suppression of J/tp mesons produced in ultrarelativistic heavy-ion 
collisions (URHICs) has long been suggested as a signature of Quark-Gluon Plasma 
(QGP) formation [[I]. However, experiments at the CERN Super Proton Syn- 
chrotron (SPS) and the Brookhaven Relativistic Heavy-Ion Collider (RHIC) show a 
very comparable suppression of J/^'s, despite much larger energy densities believed 
to be created at RHIC compared to SPS. Such a result has, in fact, been predicted 
in approaches [ [2| where suppression mechanisms are supplemented by secondary 
charmonium production, e.g., via the coalescence of c and c quarks at hadronization 
transition [ El SI- However, significant uncertainties in these approaches remain, 
most notably due to the total charm cross section and the degree of thermalization 
of the charm-quark spectra. It is therefore essential to find additional means to dis- 
criminate suppression and regeneration mechanisms. Transverse- momentum (pt) [ 
M H El [SI |5] and rapidity (y) [E] spectra, as well as elliptic flow [El d EH EC] of 
J/^'s are hoped to provide such tools. 
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In the present paper we first discuss our recent extension [ [T^] of the two- 
component model [ [51 1 1 3] to compute pt spectra, augmented by predictions for 
elliptic flow, V2{pt), in Sec. [5] In Sec. [3] we evaluate additional mechanisms relevant 
to J/ip production at high p t >5GeV/c. Sec. [4] contains our conclusions. 

2. p t Spectra and Elliptic Flow in the 2-Component Approach 

Within the two-component approach, the p t spectra of charmonia = J/i/j, Xci 
in URHICs are decomposed into contributions from direct production and regen- 
eration [ [T^] • The former is based on initial spectra in elementary N-N collisions 
including a Cronin effect (i.e., nuclear p t broadening implemented via a Gaussian 
smearing) and primordial nuclear absorption estimated from p-A collisions; the re- 
sulting spectra serve as an input to a Boltzmann transport equation to compute 
subsequent suppression in an expanding thermal fireball through the QGP, mixed 
and hadron gas (HG) phase. The charmonium dissociation rates in the QGP ac- 
count for in-medium reduced binding energies and are therefore calculated using 
inelastic quasifree cross sections [ [2] (rather than gluo-dissociation [ [14]), while 
HG dissociation is estimated using SU(Nf=4) effective theory. We also take into 
account the "leakage" effect [ US] [16], i.e., charmonia traveling outside the fire- 
ball boundary before freeze-out are not subject to dissociation. The leakage effect 
reduces suppression primarily for high-p t charmonia. The inclusive yield of the re- 
generation component is computed from c-c coalescence using detailed balance (i.e., 
the inverse of quasifree dissociation) within a momentum-independent thermal rate 
equation [ 13]. Regeneration processes are suppressed by a schematic relaxation- 
time factor to simulate incomplete kinetic equilibration of the charm quarks. The p t 
spectra of regenerated ^P's are approximated by a local thermal blastwave distribu- 
tion [ [17] at the hadronization transition based on the same fireball as used for the 
direct component. The input charm cross section, cr ec =570pb, is consistent with 
recent PHENIX data [ [18] . Our approach yields a reasonable overall description of 
the nuclear modification factor, RAA{N par t), for inclusive J/tp production, as well 
as its p t dependence in Pb-Pb and Au-Au collisions at SPS and RHIC, respectively [ 
IT2"] . At SPS and for peripheral collisions at RHIC, the regeneration contribution is 
small, but it becomes significant for p t <2-3 GeV for semi-/central collisions Au-Au 
collisions at RHIC, see, e.g., the left panel of Fig. [1] 

We now apply our approach to estimate the elliptic flow of J ftp's according to 

v 2 ( Pt ) = f dir vf r { Pt ) + r 9 v r 2 e9 {pt) , (i) 

where the fractions f dlr and f reg follow from the decomposition provided by the 
two-component model for 20-40% central Au-Au at RHIC. For the underlying v 2 
spectra we take guidance from the suppression calculations in Refs. [H] [TUj for the 
direct component, (pt), and from the coalescence models in Refs. [El [H] for the 
regenerated component, cf. the right panel in Fig. [I] The weighted sum is shown by 
the green band reflecting the uncertainties in the input v 2 of the 2 components. It is 
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Fig. 1. (Color online). Left panel: 2-component approach for J /tp Raa{pi) in 
20-40% Au-Au at RHIC, compared to PHENIX data [115]. Right panel: estimated 
V2(pt) based on the partition in the left panel. The band represents uncertainties 
due to differences between two independent input V2(pt) calculations for each of the 
two components. 



somewhat sobering to find that the resulting V2{pt) is small, around 2-3%, over the 
entire pt range. The reason is that a large V2 for the coalescence component only 
develops above pt—2 GeV, where, however, its weight is already small. On the other 
hand, in the low p t region, where the coalescence component carries a significant 
weight, its i>2 is still small. The V2 of the direct component, which solely derives from 
the path length dependence relative to the reaction plane, is always small. This puts 
a rather high requirement on the accuracy of future measurements of J /tp elliptic 
flow in order to extract information about J jip production mechanisms. 

3. Charmonium Production at High p t 

The predictions of our approach for p t spectra in 200 AGeV Cu-Cu collisions are 
compared to data in the left panel of Fig. [2] (solid lines). While the agreement is 
fair up to p t — 4 GeV, at higher momenta a significant discrepancy develops. In the 
following, we will investigate several potential sources of this discrepancy related to 
high-pi production which have not been included in our baseline calculations [I12j. 

First we check whether the increasing trend of R,AA(pt) can be attributed to a 
stronger Cronin effect, which at RHIC energies is currently not well constrained 
by d-Au data, cf. the right panel of Fig. [5] In our baseline calculations, we 
used the same broadening as at SPS, according to (Pt)AA = (Pt)pp + a gNL with 
a g Ar=0.1GeV 2 /fm (L is an average gluon path length prior to the hard process for 
charm production). With a larger a 9 jv=0.4GeV 2 /fm, which might still be accept- 
able for the d-Au data, the high p t regime in the Cu-Cu spectra indeed improves 
(dashed lines in the right panel of Fig. [2]), which is, however, somewhat at the ex- 
pense of the low-p t region. Clearly, more accurate d-Au data will be very valuable 



4 



X.Zhao and R.Rapp 



STAR Cu+Cu 0-W!' i .'STAR :+p 



1/2 



< 

DC 



\{\ O PHENIX Cu+Cu 0-6fWPHENIXp+p 

! # PHENIX Cu+Cu 0-60%/STARp+p 

- □ PHENIX Cu+Cu 0-20%/PHENlX p+p 
" ■ PHENIX Cu+Cu 0-20%/STAR p+p 

" 2-ComponL'ur Approach Cu+Qi 0-6OS a_gN=0.1GcV"2/t'm 

" 2-ComponeHl Approach Cu+Cu 0-20% a_gN=0.1GeV A 2/t'm 

- ■ - ■ 2-ComponL'ur Approach Cu+Cu 0-60% a_gN=0.4GcV"2/t'm 
■ 2-CoiiipoiiL'ui Approach Cu+Cu 0-207c a_gN=04GeV ft 2/fro 



200 A GeV _ 




d+Au s ~-200AGeV Minibias □ PHENIX 



2 4 6 8 

Transverse Momentum p T (GeV/c) 



4 6 
P. (GeV) 



10 



Fig. 2. (Color online). Left panel: J/ip RAA{pt) in Cu-Cu at RHIC within our 
2-component approach using a g jv=0.1,0.4 GeV 2 /fm (solid, dot-dashed lines), com- 
pared to PHENIX [[20] and STAR data [ 21j. Right panel: Cronin effect for J/tp 
RdA(pt) implemented via Gaussian smearing with a 5 Ar=0.1,0.4 GeV 2 /fm (solid, 
dashed line), compared to PHENIX d-Au data [[22]. 



to constrain the initial nuclear p t broadening. 

Second, we consider the following two effects primarily relevant at high p t : (1) 
Finite formation times [[TBI Ell [24] and (2) Bottom feeddown. Concerning (1), one 
expects reduced dissociation cross sections for a "pre-hadronic" cc pair relative to a 
fully formed charmonium due to a finite formation time, Tf , required to build up the 
hadronic wave function. For a schematic estimate we parameterize the evolution of 
the pre-hadronic dissociation rate as 

r c c-r*r/r} ab , t <r l f ab = T f m t /m* (2) 

with r*: (nuclear, partonic or hadronic) dissociation rates for a formed charmo- 
nium, r: fireball proper time, 17=0. 89(2. 01, 1.50) fm/c: formation time of J/ip(Xc, 
■!//) in its rest frame, and nit = +Pt ) 1 ^ 2 - In essence, (pre-) charmonium dis- 
sociation rates acquire an additional momentum dependence through Lorentz time 
dilation, being reduced at high p t . Also note that the longer formation times of 
Xc and 1// imply less suppression relative to J/tp, quite contrary to standard disso- 
ciation and statistical regeneration mechanisms, thus constituting a rather unique 
signature of the formation time effect. Concerning (2), the left panel of Fig. [3] 
shows recent data on the B — > J /ip feeddown fraction in elementary p-p(p) colli- 
sions, which is quite significant. As an estimate of this contribution, we use the 
Tevatron data [ [25] with an extra uncertainty of ±50% and replace that part of 
our direct component. Combining formation time effects and _B-meson feeddown 
we obtain as our current best estimate the green band in the right panel of Fig. [3] 
which is not far off the experimental data. 

Finally, we compare the impact of different QGP dissociation mechanisms on 
J/ijj pt spectra. In addition to quasifree destruction, we have employed the tra- 
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Fig. 3. (Color online). Left panel: Tevatron [[25] and STAR data for the ratio 
of J/ip from B-meson feeddown to inclusive J/ip [US]. Right panel: J/ip R,AA(pt) 
within the 2-componcnt approach including formation time effects and .B-meson 
feeddown (green band) in 0-60% Cu-Cu, compared to our previous results (solid 
lines) and RHIC data [[HUH]. 



ditional gluo-dissociation [ [Tl] (using vacuum binding energies) , which exhibits a 
significantly softer 3-momentum dependence, cf. left panel of Fig. [|] This indeed 
leads to somewhat harder p t spectra (see right panel of Fig. QJ, but we note that 
the use of vacuum J/tp binding energies in the QGP may not be very realistic. 
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Fig. 4. (Color online). Left panel: Comparison [[12] of the momentum dependence 
of quasifree (solid lines) and gluo-dissociation rates (dashed lines) for J/ip at differ- 
ent temperatures. Right panel: J/ip Raa{pi) with gluo-dissociation (dashed lines) 
compared to quasifree destruction (solid lines) and RHIC data [[20] [21]. Formation 
time effects and B meson feeddown are not included. 
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4. Conclusions 

Employing a 2-component model extended to calculate charmonium p t spectra in 
URHICs, we have computed J ftp elliptic flow and studied additional effects relevant 
at high p t . We find a rather small V2 due to the smallness of the regeneration 
component at high p t . A fair description of available high-p t spectra seems to be 
possible upon inclusion of formation time effects and J/ip's from B-meson feeddown. 
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